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Summary. — We report on a recent calculation of single-inclusive high-pr jet 
production in unpolarized and longitudinally polarized pp collisions at RHIC, in- 
vestigating the effect of the algorithm adopted to define the jets on the numerical 
results for cross sections and spin asymmetries. 
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1. Introduction 

Jets arc important tools in QCD for investigating the partonic substructure of hadrons 
and interactions among partons. There is no unique way to define a jet. It is thus 
important to compare and contrast the different available algorithms. The jet algorithms 
can be divided into two broad classes; (i) successive combination [1]: in this scheme, one 
defines a distance between a pair of objects and a beam distance for every object as 
follows: 



(1) 



d i:j = min(fc^ 



p2 



k 2p 



dij is called the distance between two particles i and j and diB is the distance between 



the beam and the particle; k t .i is the transverse momentum of the 
respect to the beam direction and 



-th particle with 



(2) 



At each step, the smallest of all distances is determined. If it is a beam distance, the 
object is called a jet and is removed from the event; otherwise the two objects j, k are 
combined into a single one. Examples of successive combination algorithms are the kt 
algorithm [2], where p = 1, and the anti-fc t algorithm [3] for which p = — 1. 
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(ii) Cone algorithms: in these algorithms the jet is defined in terms of stable cones 
as circles of fixed radius in the rj-<p plane, such that the sum of the 4-momenta of the 
particles in it points in the direction of the center of the cone. One defines the jet by all 
particles j that satisfy [4] 

(3) Rlj ^ ( VJ - V] ) 2 + (cjij - h) 2 <R 2 , 

where r\j and 4>j are the pseudo-rapidity and azimuthal angle of the jet, respectively. 
Higher order QCD corrections are important as the dependence on the factorization and 
renormalization scales is expected to be reduced when the corrections are included. In 
the case of jet production, higher order corrections are particularly important, as only 
at NLO the QCD structure of the jet starts to play a role in the theoretical description 
of the process. In fact, some of the popular cone algorithms are known to be collincar 
and infra-red unsafe at NNLO or when multiple jets are considered. 

Single inclusive large pt jets in longitudinally polarized pp collisions at RHIC are 
important tools to gain access to the polarized gluon distribution in the nucleon. The 
cross section for single inclusive jet production at RHIC has been calculated at NLO using 
a Monte-Carlo technique [5] in the cone algorithm. However a largely analytic technique 
was developed in [6] in the limit when the cone opening is relatively small (small cone 
approximation). This is advantageous because it leads to much faster and more efficient 
computer codes as the singularities in the intermediate steps cancel analytically and 
one does not have to treat them through delicate numerical techniques. The basis of 
such an analytic calculation is the observation that the inclusive jet production proceeds 
through the same partonic subprocesses as single inclusive hadron production and it 
is possible to convert an NLO cross section for single inclusive hadron production to 
the one for jet production. The main difference between the two cases is the fact that 
in single inclusive hadron production, one integrates over the full phase space of the 
unobserved partons. This leads to collincar singularities, which are absorbed in the 
parton to hadron fragmentation functions. In contrast, for a jet, final-state particles 
that move in roughly the same direction will jointly produce the jet. This makes the 
cross section more inclusive, and (for a proper jet definition) final state singularities 
must cancel. The cross section for single inclusive hadron production can however be 
transformed into that for single inclusive jet production [6]. In the limit of small cone 
size, this transformation can even be performed analytically. 

We have recently extended the above analytic technique to the more widely used 
successive combination schemes (for example, k t or anti-/c t ), assuming again that the jet 
parameter R used to define the distance between two objects in this algorithm is not too 
large [7]. When systematically expanded around R = 0, the dependence of the partonic 
cross sections on R is of the form _41ogi? + B + 0(R 2 ). The coefficients A and B are 
calculated analytically, and the remaining terms 0(R 2 ) and beyond are neglected. Wc 
refer to this approximation as "Narrow Jet Approximation" (N J A) . The N J A gives a very 
accurate description of the single inclusive jet cross sections at RHIC, Tevatron and even 
at the LHC [7]. It turned out that the cross sections for single inclusive jet production 
in the cone and the successive recombination algorithms differ by calculable finite terms. 
We have also given numerical estimates of the cross section both for unpolarizcd and 
longitudinally polarized collisions at RHIC, and examined the effect of the choice of jet 
algorithm on the double longitudinal spin asymmetry. Here wc give a brief report of [7] . 



JET PRODUCTION IN PP COLLISIONS: DEPENDENCE ON JET ALGORITHM 



3 



2. Cross section for single inclusive jet production in pp collisions 

We consider single-inclusive jet production in hadronic collisions, pp jet JT, where 
the jet has a transverse momentum ptj, rapidity r]j, and azimuthal angle <pj. Note 
that on top of the choice of jet algorithm one also has to define how objects are to be 
merged to form the jet. We choose to define the four-momentum of the jet as the sum 
of four- momenta of the partons that form the jet for both algorithms ( "E recombination 
scheme" [4]). 

In order to calculate the single-inclusive jet cross section at NLO, we start from the 
NLO single-parton inclusive cross sections dcr a b^, c x , relevant for single- inclusive hadron 
production process pp — > hX and analytically known. For a jet cross section, the observed 
final state should not be given by parton c only, but by partons c and d jointly, when 
the two are close to each other (as two partons together can form the jet). In order to 
calculate this one first considers a "jet cone" characterized by a jet parameter R around 
the observed parton c and notices that in the NLO single-parton inclusive cross section 
there is a configuration where an additional parton d is inside the cone (we use the term 
"cone" for simplicity, the considerations apply to any jet definition). One subtracts these 
contributions and replaces them by terms for which partons c and d are both inside the 
cone and form the observed jet together. For a given partonic process ab — > cde we then 
have, 



Here dbj is the single-parton inclusive cross section where parton j is observed (which 
also includes the virtual corrections), d<7j(k) is the cross section where parton j is observed 
but parton k is also in the cone, and d&jk is the cross section when both partons j and 
k are inside the cone and jointly form the jet. One has to note that the single-parton 
inclusive cross section contains a subtraction of final-state collinear singularities in the 
modified minimal subtraction (MS) scheme. One has to perform an MS subtraction also 
of the singularities in the rfoj(fc) + d&k(j) — d&jk- 

The difference between the cone and k t type algorithms resides entirely in the d&jk, 
which can be calculated analytically in the NJA. The reason for this difference is as 
follows: for the fct-type algorithms the two partons j, k are merged into one jet if their 
distance defined in (1) is smaller than their respective beam distances d\B and djB defined 
in (1). For dbjk this has to hold, and we arrive at the condition 

(5) Rj k < R 2 for k t — type algorithms, 

with Rjk defined in Eq. (2). This condition is true for all fcj-type algorithms. Whereas 
in cone algorithm, Eq. (3) is valid: 



(4) 



d<J a b-> jctX 



= [da c — da c ( d ) — da c ( e )\ 
+ [d&d — da d ( c ) — d<7d( e )] 
+ [da e - da e{c) - da e{d) ] 
+ da c d + da ce + dbde- 



(6) 



R]j < R 2 A Rlj < R 2 for cone algorithm. 



We find that the difference in the cross sections calculated for the two algorithms is finite, 
as it must be. 
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Fig. 1. - The ratio ^ a i go at RHIC for VS = 200 GeV (left) and VS = 500 GeV (right), for the 
spin-dependent case. Results are shown for two different values of the jet parameter R. We have 
chosen the factorization and renormalization scales as hf = fin =Pt,- 



3. Numerical Results 

Next, we present some numerical results for single-inclusive jet production cross sec- 
tions and spin asymmetries in pp collisions at RHIC. We use the CTEQ6.6M parton 
distributions [8] for the unpolarized cross section and the "DSSV" helicity parton distri- 
butions of Ref. [9] for the polarized case. We define the ratio 

[) alSO " mA)a/d PTj d VJ ] conc ' 

where the jet parameter R is the same for both cross sections. 

Figure 1 shows the ratio 7?. a igo for polarized collisions at RHIC, calculated for the 
factorization and renormalization scales as [ip = /in = pt j1 as a function of pxj in bins 
of ptj , for R = 0.4 and 0.7. We present results for two values of c. m. s. energies at 
RHIC, \/S = 200 GeV (left) and \/S = 500 GeV (right) . TZ a i go is around 90% at high 
PT, , but deviates largely from one in the bin around ptj = 12.5 GeV. The reason for this 
is that for the DSSV set of parton distributions the polarized jet cross section changes 
sign around pxj = 10 GeV. Depending on the jet algorithm, the zero is at a slightly 
different value of px- This shows that in regions where the polarized cross section is very 
small it is very sensitive to the choice of jet algorithm. 

The double longitudinal spin asymmetries All at RHIC arc defined by 

(8) A LL = d2A<J / d PTj dr )J 

d 2 <7 j 'dpTjdrjj 

For the denominator we use the spin-averaged cross sections and for the numerator the 
polarized ones. The results are shown in Fig. 2. As one can see, the asymmetries are 
quite insensitive to the jet algorithm chosen, and also to the value of the jet parameter 
R for all values of pT r where the asymmetry is sizable. Our results are useful for the 
analysis of the data on the double longitudinal spin asymmetry in single inclusive jet 
production by the STAR collaboration at RHIC [10]. 
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Fig. 2. - Double-longitudinal spin asymmetries All at RHIC, for s/Ti = 200 GeV and \/~S = 500 GeV 
and various jet definitions. We have averaged over \r]j\ < 1. The factorization and renormalization 
scales have been chosen to be hf = /-in — Ptj- 
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